A Monte-Carlo and Self-Consistent Field calculations of encapsulated spherical 
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We present the results of extensive numerical self-consistent field (SCF) and 3-dimensional off- 
lattice Monte Carlo (MC) studies of a spherical brush confined into a spherical cavity. The monomer 
density profile and the cavity pressure have been measured in systems where curvature of the cavity 
has an important effect on the polymer brush conformation. A direct comparison between the SCF 
and MC methods reveals the SCF calculation to be a valuable alternative to MC simulations in 
the case of free and softly compressed brushes. In the case of strongly compressed systems we have 
proposed an extension of the Flory theory for polymer solutions, whose predictions are found to 
be in good agreement with the MC simulations and has the advantage of being computationally 
inexpensive. In the range of high compressions, we have found the monomer volume fraction v to 
follow a scale relationship with the cavity pressure P, P '-^ v" . SCF calculations give a = 2.15±0.05, 
close to des Cloiseaux law (a — 9/4), whereas MC simulations lead to a = 2.73 ±0.04. We conclude 
that the higher value of a obtained with MC comes from the monomer density correlations not 
included in the SCF formalism. 

PACS numbers: 



o 



> 

o 
o 
o 

-(— • 



I 

o 



% 



I. INTRODUCTION 

The structure and dynamics of polymer chains ter- 
minally anchored or end-grafted to a surface is a very 
interesting problem from both, theoretical and experi- 
mental points of view, due to the inherent complexity of 
having chain molecules located in constrained geometri- 
cal environments. These polymer systems are relevant 
in many areas of polymer science and technology [ll Q], 
such as chromatography, surfactants, lubrication, adhe- 
sion, stabilization of colloidal systems, or even as po- 
tential drug carriers[3, H, IM S' Properties related to 
these systems are also of interest in the synthesis of 
nanoparticlesfTj; spherical brushes formed by colloidal 
poly- (methyl methacrylate) spheres with a grafted layer 
of poly(12-hydroxy stearic acid) (PMMA-PHSA) have 
served as a model to emulate hard-spheres minimizing 
the effect of van der Waals forces [8| • More recently, the 
use of the light-scattering properties of coated spherical 
particles has been found to be a useful non-destructive 
way to probe and size systems ranging from blood cells 
to paper whiteners Q- 

Former theoretical studies by de Gennes and 
Alexander [lOl 111 . Semenov [ll, Milner-Witten- 
Gates [13|, [1^ and Zhulina-Priamitsyn-Borisov 15J fo- 
cused on polymer chains grafted onto planar surfaces, 
and their predictions have been well supported by Monte 
Garlo[l^[l3illSll3 and molecular dynamics ^0!| simula- 
tions. 

The introduction of curved interfaces result in poly- 
mer brush structures whose properties differ significantly 
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from those expected from flat interfaces. Whereas the 
problem of concave surfaces, in which the surface curves 
towards the polymer chain, was solved by Semenov |2]]|. 
for convex surfaces it has not been possible to find a po- 
tential form that provides a self-consistent solution lead- 
ing to physical solutions. In this case, the increased vol- 
ume available to the stretched polymer, as it moves away 
from the interface, is responsible of a rich physical be- 
havior. For instance, colloidal silica spheres with grafted 
alkane chains undergo a sol-gel transition when dispersed 
in hexadecane solution|2^. 

On the other hand, experimental studies on the rhe- 
ological properties of spherical brushes appear rather 
scarce, certainly due to the str ong difficulties in the con- 
trol of the grafting procedure l23[. Nonetheless, several 
studies using SANS 24,'25], ESR 26, "2?] or NMR [Hill 
technics have dealt with the physical chain-properties of 
terminally attached homopolymers onto colloidal parti- 
cles. 

^From the theoretical point of view, there have 
been several attempts to determine the interaction be- 
tween spherical brushes. Daoud-Gotton^C^ and Witten- 
Pincus 31j in the limiting case of star shaped polymers; 
Borukhov and Leibler|33| using the Derjaguin approxi- 
mation valid for small curvature effects and short poly- 
mer chains. Several numerical studies have been de- 
voted to the study the properties of single non- interacting 
brushes [3^ 133. 1351 and the interaction between two 
spherical brushes [Sa, [23 ■ 

Whereas the previous studies have analysed the prop- 
erties of a single or two interacting spherical brushes, 
the understanding of the more complex and interesting 
problem of a system of colloidal spherical brushes in so- 
lution remains untouched. Particularly, in the limiting 
case of high densities, the colloidal brushes are subject 
to an isotropic pressure. In such case, the interaction of 



a colloidal particle with the rest of the system can be 
modelled, in a first approximation, by a single spherical 
brush confined inside a spherical cavity. This model has 
the advantage of reducing the high computer cost. On 
the other hand, this model turns out to be relevant to 
the study of the properties of encapsulated dendrimers, 
liposomes, vesicles containing nanoparticles with grafted 
chains, and might be of relevance in the synthesis of 
polymer-grafted metal nanoclusters inside small material 
cavities or molecular cages. 

The purpose of the present paper is to present the re- 
sults of extensive MC and SCF simulations of a single 
colloidal brush confined inside a spherical cavity wall of 
variable radius. Polymer chains have an extent of the 
same order of the core size of the colloidal particle, such 
that curvature effects are important. We have measured 
the monomer density profile and the cavity pressure at 
different cavity radius. The results using MC and SCF 
methods have been compared at low and high compres- 
sion regimes, discussing the advantages of each method. 

The rest of the paper is organized as follows: in section 
II we describe the numerical procedures used to compute 
the monomer density profile and the cavity pressure; in 
section III we present a detailed analysis of the numerical 
calculations, this section also includes an extension of the 
Flory theory for a free polymer solution, as an alternative 
way to obtain the pressure inside the cavity; section IV 
concludes with a brief summary and a discussion of the 
results. 



II. NUMERICAL MODELS 
A. Monte Carlo method 

In order to simulate the interaction between a spher- 
ical brush confined inside a spherical cavity wall, we 
have used a 3-dimensional off-lattice Monte Carlo method 
(MC). We have generated the brush by homogeneously 
distributing / polymer chains grafted onto an impene- 
trable spherical surface of radius r^. The cavity wall is 
also impenetrable with a variable radius R. The polymer 
chain is represented by the pearl necklace model [33 con- 
taining TV beads of diameter a. The distance between two 
consecutive beads in the chain is set to l.lcr. In all the 
simulations we have set a = 1. The initial configuration 
of the self-avoiding polymer is randomly generated being 
the first monomer permanently anchored to the surface 
(it is never allowed to move). An schematic representa- 
tion of the system is shown in Figure ^ 

Monomers interact through a steric hard-core potential 
of the form: 



where T^ is a hard sphere potential: 
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Different polymer configurations are generated by 
changing the position of a randomly selected monomer. If 
the monomer is located inside the chain (between the first 
monomer, permanently anchored, and the last one) its 
position changes by rotating an arbitrary angle between 
and 27r around the axis connecting the previous and 
following monomers in the chain. Chain ends just per- 
form random wiggling motions. The proposed motion is 
accepted if the excluded volume interaction is preserved 
(Eq. 121). A link-cell method [33 has been implemented 
in the algorithm to efhciently check all possible monomer 
overlaps. 

Initially, the radius of the cavity wall is set to be larger 
than the usual extent of the brush to ensure that dur- 
ing the equilibration process no interaction between the 
brush and the cavity wall occurs. We define one Monte 
Carlo Step (MCS) as N x f trials to perform monomer 
moves. The spherical brush has been equilibrated typi- 
cally during 5 x 10^ MCS. After this initial equilibration 
time, magnitudes of interest are recorded every 10 MCS. 
In what follows, the expression for the force and inter- 
acting potential are given in units of ksT. 

The force that the system exerts to avoid compression 
is computed as the change in the free energy J- due to an 
infinitesimal change in the cavity radius R is given by: 
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where Z{R) = CKp(—J-/kBT) is the partition function 
of the system. Due to the hard-core structure of the 
potential, see Eq.Q, the partition function 
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is equal to the volume Vt{R) of all possible polymer chains 
configurations compatible with the steric requirements. 

Let $7(7 (i?) C Vt{R) be the subset of configurations 
compatible with a reduction of the radius cavity by an 
amount 6R. Similarly, let VIe{R - 5R) C n{R - SR) be 
the set of polymer configurations oi n{R — SR) in which 
it is possible to increase the radius of the cavity by an 
amount (5 i?. The compression probability is just Pc(i?) = 
^ciR)/^iR)i and the probability that an expansion can 
be done is PEir - SR) = 17i=;(i? - SR)/fl{R - SR). Since 
no steric requirement can prevent the expansion of the 
cavity, Pe{R) = 1 for all values of i?. The subsets il,c{R) 
and flsiR — SR) are in one to one correspondence, con- 
sequently 
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Because of the finite step 5R used in computer simu- 
lations, we approximate 
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The compression probabihties are obtained by span- 
ning directly this probability for a given cavity size R. 
Then, we force to compress the system to a new cavity 
size R — SR. The system is equilibrated again during 10^ 
MCS before new measurements are taken. 

The monomer radial density (f> is defined as usual: 



Anr'^dr' 
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being i>(r) the number of monomers located within a dis- 
tance between r and r + dr from the centre of the sphere. 
The definition of (j){r) is such that the following normal- 
ization holds: 



Nf 



(8) 



In a similar way, we also define the chain-ends monomer 
density, e(r), where iy{r) takes only into account the num- 
ber of free chain ends. 



B. Self Consistent Field method 

In order to compute the probability density func- 
tion (pdf) for polymer systems, it is customary to deal 
with a Schrodinger-like equation for the pdf for a single 
chain J4flj. The use of this formalism is a valid approxi- 
mation for spatial scales much larger than the polymer 
blob size. However, for encapsulated polymer brushes, 
and mainly at moderate and high compressions, the spa- 
tial scale of the cavity is comparable to the blob size. 
Thus, the use of the precedent formalism to compute the 
pdf becomes inadequate. Instead, we compute the pdf of 
a polymer chain using directly the recurrence law for the 
pdf from which the Schrodinger-like equation derives |4flj. 
We must be aware that the use of the recurrence law still 
implies some approximations. It is assumed that on each 
monomer of the polymer chain acts a potential that only 
depends on the position of the monomer in the system, 
U = U{r), therefore, bond correlations are not taken 
into account. It is also assumed that the potential is a 
function of the local monomer density, U{(j){r)), ignoring 
particle density correlations. Finally, the properties of 
the whole ensemble of chains are deduced from the pdf 
of a single chain. 

Under the above assumptions, the spherical cavity is 
discretizcd in concentrical shells of thickness dr and all 
monomers inside a shell are assumed to be equivalent. A 
polymer chain composed by N monomers is represented 
as a path of A'^ segments of length a. Each segment is 
labeled by an index r associated to the spherical shell at 
which it belongs. The pdf associated to all possible paths 



composed by n segments, being the first segment inside 
the shell h' and the last one inside the shell h, is defined 



as 
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in which ^ stands for a sum over all the hypothetical n- 
paths that join the shells h' and h. This function verifies 
G'„(h',h) = G'„(h, h'). Therefore, the pdf associated to 
a path of n + 1 segments may be written as 



„+i(h',h) = Y.^ (e-E?^i^(-0 e-^(-.+i))(io) 
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Assuming U{Tn+i) to be independent of rest of the seg- 
ments in the chain and of the starting point of the se- 
quence, the precedent equation reads: 

G„+i(h',h) = f^^G„(h',h")j e-^C^) (11) 

where ^ ^ implies a sum over all the h" shells from 
which we can get into shell h using a single segment, 
therefore, shells h" and h are at a relative distance less or 
equal to a. The above equation stands for the recurrence 
law needed to calculate the chain pdf once the potential 
U{h) is given. 

We have set the interaction potential U(i) to be pro- 
portional to the monomer concentration in shell i, (fiii): 

U{i) = Lu (j){i) (12) 

where w is the excluded volume parameter defined as|4H: 
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where V{r) is the interaction potential introduced in 
Equation Q, that for a value of cr = 1, as in the MC 
method, we obtain oj = Att/S « 4.2 

The monomer concentration 0(i) and the free-end- 
chain concentration e(i) are defined as: 
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where / is the number of chains and v^ and i^ are the 
shell indexes with radius equal to the brush core sur- 
face and the cavity wall respectively. The factor e*^*^'^ in 
the monomer concentration is introduced to avoid dou- 
ble counting of the interaction term e"*^^*^ at shell i that 
comes from splitting the pdf into two terms, one that 
ends at shell i and the other starting at the same shell. 



The change in the free energy when we reduce the cav- 
ity size, AT^ifi), is given by 



ATiii 



In 
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being 
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We have set the reference state as a system with a cavity 
wah at a distance far enough so that no chain can reach 
this wall. Therefore, the force to compress the brush, 
given a cavity size i?, is: 



F{R) = 



dR 



(18) 



Given R, N and f, an iterative process is used to ob- 
tain the pdf and the density profiles. We have iterated 
the process until self-consistency is reached. We set the 
condition for self-consistency such that the sum of the 
square differences in the density profiles coming from two 
consecutive iterative steps is less than 10~*. 

The SCF method has the great advantage of being 
three to four orders of magnitude less expensive in com- 
puter time than the MC method. In the next section we 
will show that SCF and MC calculations give similar re- 
sults for the density profiles (/)(r) and e{r) in the case of 
free polymer brushes. However, SCF predictions for the 
cavity pressure worsen for highly compressed systems. 



III. RESULTS AND DISCUSSION 
A. Density profiles 

We have performed extensive numerical calculations 
for free and encapsulated spherical polymer brushes in 
order to computed the monomer density profiles and 
the compression forces for different sets of parameters 
(i?, N^ /). The core radius of the colloidal particle where 
polymers are grafted is taken to be Tc = 5cr and is kept 
constant through all the simulations. The diameter of 
the monomers is set to a = 1. Wc have taken polymer 
chain lengths in the range of A^ = 30 to A^ = 70, and we 
have varied the number of grafted chains from / = 5 to 
/ = 75. The range of parameters (iV, /) has been chosen 
in order to obtain chain extents roughly of the same or- 
der than the diameter of the core where curvature effects 
are important. For the SCF method we have used a shell 
thickness dr = O.lcr. 

We have first studied the monomer density profiles 4>{r) 
for an unconstrained spherical brush using our MC sim- 
ulations and the SCF formalism. In Figure |21 we present 
the results for two different values of the chain length N 
and number of grafted chains / that correspond to repre- 
sentative values for soft and densely packed brushes. We 
can observe a good agreement between the MC and SCF 



calculations. The density oscillations observed at small 
r, close to the core of the colloidal particle, are origi- 
nated due to wall-effects of the impenetrable core. For 
a free spherical polymer brush Cariagno and Szleiferpa| 
computed the monomer density profile derived from a 
single-chain mean field theory. A comparison between 
the Cariagno and Szleifer data (CS) and our results is 
also included in Figure |21 The better agreement of the 
CS predictions with our MC simulations, in contrast to 
the SCF calculations, can be understood in the sense that 
the Cariagno and Szleifer formalism requires a represen- 
tative sample of chain configurations as an input data to 
solve the equations that we have generated using our MC 
method. 

In Figure |3| we compare the monomer and chain- 
end density profiles for unconstrained polymer brushes 
obtained from our MC and SCF calculations. Figure 
ISJa) stands for A^=30 and /=25, whereas Figure Elb) 
and Figure |2{c) show the results for (iV=30,/=75) and 
(iV=50,/=75) respectively. In all cases the profiles are 
roughly similar. In particular, both methods agree very 
well in the case of the chain-end density profile e(r) for 
short chains in a densely packed brush (see inset of Fig- 
ure|3Jb)). This is mainly due to the fact that the chains 
are forced to be mostly fully stretched out and den- 
sity correlations, not present in the SCF formalism, are 
not relevant. On the other hand, the results for the 
monomer density profiles 4'{r) show systematic differ- 
ences, although small, between the MC and SCF cal- 
culations. Close to the core of the colloidal brush, SCF 
results display density profiles slightly smaller than those 
obtained via MC simulations. And vice versa, in an inter- 
mediate region, the SCF method gives densities slightly 
larger than in the MC simulations. The same system- 
atic behaviour was found by Cosgrove et al.|4J|, when 
comparing MC and SCF density profiles for flat brushes. 
Cosgrove attributed these differences to the fact that MC 
simulations accounts explicitly for the excluded volume 
effect, whereas SCF accounts only approximately for this 
effect. 



B. Cavity pressure and force profiles 

We have measured the force proflle exerted by an en- 
capsulated spherical polymer brush onto the external 
cavity wall through the evaluation of the changes in the 
free energy due to an infinitesimal change in the radius of 
the cavity. Within the MC simulations, the force can be 
calculated by directly measuring the compression proba- 
bilities (Equation O; whereas in the SCF approach, once 
we have reached self-consistency, we use the pdf in Equa- 
tions [H to UHI 

Alternatively, we propose another way to compute the 
force profile and the cavity pressure that is derived from 
the Flory theory for polymer solutions 4i|| . This method 
has the advantage of being computationally inexpensive, 
and their predictions will be compared with the MC and 



SCF calculations. 



An extension of the Flory theory for polymer solutions 

In the Flory theory the osmotic pressure in a polymer 
solution can be written as: 

IV{R) ^ ^(\n{l-v) + {l-^)v + xv^^ (19) 

where V is the molar volume of the solvent, v is the vol- 
ume fraction of solute, N is the degree of polymerization 
or chain length, and x is the Flory parameter. We set 
X = which is the condition of a dry-brush. Under this 
assumption, the contributions to the free energy come 
only from the entropy associated to all possible configu- 
rations of the system. We suppose that the volume of the 
solvent and the volume fraction of the solute are given 
respectively by: 



V ^ {T-aNf) 
aNf 



(20) 
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a is the volume of a single monomer, r is the total avail- 
able cavity volume between the inner wall, represented by 
the core of the colloidal particle where chains are grafted 
and the cavity wall, thus: 



47r , 
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(22) 



At variance with the original Flory treatment, the mo- 
lar volume of the solvent V refers to the remaining space 
in the system once we have subtracted the volume occu- 
pied by the monomers, thus, it is no longer a constant 
value. 

The force to compress the cavity will be proportional 
to the area of the cavity wall and to the change in the 
osmotic pressure, thus: 



F{R) -47ri?2An(i?) 



(23) 



In Figures ^ to El we present in log- log plots the force 
profile F{R) vs. the cavity size R computed for differ- 
ent values of the polymer chain length N and number of 
grafted chains /. In each figure we include the results 
coming from the MC simulations (circles), the predic- 
tions of the SCF theory (dashed lines) and results de- 
rived from the application of the extended Flory theory 
(crosses). Figure0]concentrates on the results derived for 
short polymer chains {N = 30); Figure|Slfor intermediate 
chain lengths {N — 50); and Figure for long polymers 
{N = 70). In all the cases studied we have used the same 
fitting constant to adjust the predictions of the extended 
Flory theory (see Eq. I23|l . and we have taken, as a ref- 
erence state, a cavity size R* at which \n{F{R*)) -^ — cx) 
in the MC simulations. 



A direct comparison between SCF and MC force pro- 
files shows a rather good agreement for weakly com- 
pressed systems. However, systematic differences are ob- 
served for intermediate and high compression values. In 
the intermediate region we found the SCF forces to be 
larger than the ones derived from the MC simulations, 
whereas for high compressions it is the MC force the one 
that becomes larger than the SCF outcome. 

The mismatches observed between the SCF and MC 
results are due to a twofold effect. For small cavity 
sizes or highly compressed systems, the assumption in 
the SCF model of a linear dependence of the mean field 
potential with the monomer density (see Eq. 1121) breaks 
down. In fact, under this assumption, the SCF formal- 
ism allows a cavity size smaller than the volume occupied 
by the polymers without requiring an infinite force. On 
the other hand, for intermediate compression values, the 
larger forces obtained with the SCF formalism are origi- 
nated in an overestimation of the monomer interactions. 
The SCF method does not include the effect of monomer 
correlations, thus it allows higher average densities in the 
system than the ones found in the MC simulations. As 
a consequence, stronger repulsions between the polymers 
take place and a higher force is required to compress the 
brush. 

The predictions of the extended Flory theory are 
found, remarkably, despite its simplicity, to be in a very 
good agreement with the results of the MC simulations 
in the intermediate and high compression regimes. How- 
ever, the force is overestimated for weakly compressed 
systems. This result is easily explained since the Flory 
theory was formerly developed for free polymer chains. 
As the cavity size grows, the difference between a system 
of grafted chains and a polymer solution becomes evi- 
dent and, as it is expected, the interaction of the grafted 
chains with the outer surface is much weaker than the 
one coming from a polymer solution. 

It is worth to notice that the MC data can be fitted 
remarkably well with the extended Flory theory for inter- 
mediate and high compression values, and with the SCF 
formalism for weakly compressed systems. 

We have analysed the relationship between the 
monomer volume fraction v (see Eq. I21|) and the pres- 
sure exerted on the cavity wall P defined as the force to 
compress the system divided by the area of the cavity. 
The results for the MC and SCF calculations are shown 
in log-log plots in Figured and Figure |H| respectively. In 
both cases, we find a complex behavior of the monomer 
volume fraction with the cavity pressure for weakly com- 
pressed systems that depends on the different values of 
the polymer chain length N and number of grafted chains 
/. But, for increasing values of the cavity pressure v and 
P follow a power-law of the form P r^ v" independent 
of N and /. The best fit to the numerical data gives 
a slope of a = 2.73 ± 0.04 for the MC simulations and 
a = 2.15 ± 0.05 for the SCF results. The exponent ob- 
tained for the SCF data is very close to the des Cloiseaux 
power law (a = 9/4) found in semi-dilute polymer solu- 



tions 4Cf| . For large monomer concentrations the polymer 
theory predicts that all thermodynamic properties must 
reach values that are independent of the degree of poly- 
merization, as we have observed in both MC and SCF 
methods. We must be aware that des Cloiseaux law is 
deduced from scaling arguments and neglects non-linear 
dependences on the concentration. This fact might ex- 
plain the agreement between SCF and des Cloiseaux law, 
whereas MC results follows a power law which exponent 
is larger due to monomer density correlation effects. 

Differences between force profiles derived from SCF 
formalisms, and those obtained from other methods that 
account for chain interdigitation and correlations be- 
tween the nearest-neighboring bonds are also referenced 
in several works. For instance, Ruckenstein-Li'45| (to 
be referred as RL) compared the experimental force pro- 
file of two interacting crossed cylinders bearing grafted 
polymer chains, with the numerical data obtained with a 
generator matrix formalism and from SCF methods. RL 
found the matrix formalism to provide a better agree- 
ment with the experimental data than the SCF results. 
In most cases, the force profiles derived by RL with the 
matrix formalism and the SCF method are rather simi- 
lar to the ones we found comparing the MC simulations 
and the SCF calculation. This reinforces our presump- 
tion that interdigitation and monomer correlations are 
responsible of the observed differences between the MC 
and SCF data. Additionally, the fact that mean field cal- 
culations overcounts the segment-segment interactions in 
the evaluation of the free energy, as mentioned above, has 
already been noticed by Lin and Cast jig . 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In this paper we have studied the behavior of an en- 
capsulated spherical brush inside a spherical cavity. This 
system is a first approximation to model the interactions 
in a solution of colloidal particles bearing grafted poly- 
mer chains onto their surface. This model is particu- 
larly relevant for moderate and high density values where 
the colloidal brushes are subject to an isotropic pressure, 
and might be of relevance to the behavior of encapsu- 
lated dendrimers, liposomes or vesicles containing poly- 
mer brushes as well. 

We have measured the monomer density profile and the 
cavity forces through extensive 3-dimensional off-lattice 
Monte-Carlo simulations and using a Self-Consistent 
Field formalism. In the latter case, we have used di- 
rectly the probability density function recurrence law for 
the propagator GAr(r, r'), avoiding the length scale ap- 
proximation involved in self-consistent field methods that 



uses Schrodinger-like equations. Alternatively, we have 
proposed a theoretical description based on an extension 
of the Flory theory for polymer solutions to compute the 
pressure inside the cavity. 

A comparison of the predicted forces exerted by the 
polymer brush onto the cavity surface among the differ- 
ent methods reveals the following: i) For weakly com- 
pressed systems, MC and SCF data show a rather good 
agreement. However, the force is overestimated in the 
extended Flory theory. This difference arises since the 
Flory theory was developed for free polymer chains and 
not for polymer brushes, thus, its prediction is not physi- 
cally relevant when the cavity wall is located at a distance 
larger than the typical brush extension, ii) For interme- 
diate and highly compressed systems the MC data agrees 
reasonably well with the results derived from the Flory 
theory. In the intermediate regime, it is the SCF formal- 
ism that overestimate the force. This behavior can be 
easily explained since the SCF method does not account 
for monomer correlations, allowing higher monomer den- 
sities and thus, higher forces are required to compress 
the brush. On the other hand, for highly compressed 
systems, the linear dependence of the mean field poten- 
tial with the monomer density turns out to be inadequate 
(the repulsion between monomers allows a volume reduc- 
tion of the system beyond the own excluded volume of 
the monomers at a finite energy cost) leading to lower 
forces than the ones derived from MC simulations. 

We have found a power law relationship between the 
monomer volume fraction and the cavity pressure, P ~ 
v". SCF data gives a slope of a = 2.15, very close to 
the des Cloiseaux law derived for semi-dilute polymer 
solutions. On the other hand, the MC simulations pro- 
vide a larger exponent a = 2.73 that is originated in the 
monomer correlations not present in the previous models. 

The study of a spherical brush constrained by an 
isotropic steric wall is directly related to the problem of 
a dense solution of colloidal particles with grafted chains 
onto their surface where chain interdigitation among 
brushes is not favoured. We expect these results to stimu- 
late further theoretical and experimental studies towards 
the understanding of the behaviour of colloidal particle 
systems in solution. In this sense, the computed force can 
be directly related to the measure of disjoint pressures by 
using the expressions derived by Evans and Napper|47|. 



Acknowledgments 

Financial support from the Spanish MCyT grants nos. 
BFM2000-1108 and BFM2001-0341-C02-01 are acknowl- 
edged. 



[1] Napper, D. H. Polymeric Stabilization of colloidal disper- 
sion: Academic, London, 1983. 



[2] Russel, W. B; Saville, D.A. and Schowalter, W.R. Col- 



loidal Dispersions; Cambridge University Press, Cam- 
bridge, 1989. 

[3] Hsu, W. P.; Yu R. and Matijevic, E., J. of Colloid Inter- 
face Sci. 1993, 156, 36. 

[4] Meyer, R. A., Appl. Opt. ,1979, 18, 585. 

[5] van Zanten, J. H. and Monbouquette, H. G., J. Colloid 
Interface Sci. 1991, 146, 330. 

[6] Kuhl, T. L. ; Leckband, D. E.; Lasic D. D.; and Is- 
raelachvili, J. N., Biophys. J. 1994, 66, 1479. 

[7] Mangeney, C. ; Ferrage, P.; Aujard I. et al, J. Am. Chem. 
Soc, 2002, 124, 5811. 

[8] Phan, S. ; Russcl, W. B.; Cheng, Z.; Zhu, J.; Chaikin, 
P. M.; Dunsmuir J. H.; and Ottewill, R.H., Phys. Rev. E 

1996, 54, 6633. 

[9] Quirantes, A.; and Delgado, A. V., J. Phys. D: Appl. Phys. 

1997, 30, 2123. 

[10] de Gennes, P. G., J. Phys. (Paris), 1976, 37, 1443; Macro- 
molecules, 1980, 13, 1069; C. R. Acad. Sci. (Paris), 

1985, 300, 839. 
[11] Alexander, S., J. Phys.(Paris), 1977, 38, 983. 
[12] Semenov, A. N., Sov. Phys. JETP, 1985, 61, 733. 
[13] Milner, S. T.; Witten T. A.; and Gates M. E., Macro- 
molecules 1998, 21, 2610; Europhys. Lett., 1988, 5, 413; 

Macromolecules 1989, 22, 853. 
[14] Milner, S. T.; and Witten, T. A., J. Phys. (Paris), 1998, 

49, 1951. 
[15] Zhulina, E. B.; Borisov O. V.; and Priamitsyn, V. A., J. 

Colloid Interface Sci., 1990, 137, 495. 
[16] Chakrabarti, A.; and Toral R., Macromolecules, 1990, 

23, 2016. 
[17] Lai, P.Y.; and Binder, K., J. Chem. Phys., 1991, 95, 

9288. 
[18] Toral, R.; Chakrabarti, A.; and Dickman, R., Phys. Rev. 

E, 1994, 50, 343. 
[19] Chakrabarti, A.; Nelson P.; and Toral, R., Phys. Rev. 

A, 1992, 46, 4930; J. Chem. Phys. 1994, 100, 748. 
[20] Murat, M.; and Crest, G.S., Macromolecules, 1989,22, 

4054; Phys. Rev. Lett. 1989, 63, 1074. 
[21] Semenov, A. N.; Sov. Phys. JETP, 61, 733, (1985) (Zh. 

Eksp. Teor. Fiz., 1985, 88, 1242. 
[22] Grant, M. C; Russel, W. B., Phys. Rev. E, 1993, 47, 

2606. 
[23] Nommensen, M. H.; Duits, G.; van den Ende, D.; and 

Mellema, J., Phys. Rev. E, 1999, 59, 3147. 
[24] Cosgrove, T.; Crowley, T.L; Vincent, B.; Barnett K.G.; 

and Tadros, Th. P., Faraday Discussions, 1982, 16, 101. 
[25] Beaufils, J. P.; Hennion, M.C.; and Rosset, R., J. de 

Physique,1983, 44, 497. 
[26] Hommel, H.; Legrand, A.P; Balard, H.; and Papirer, E., 

Polymer, 1983, 24, 959. 
[27] Hommel, H.; Legrand, A. P., J. le Courtier and J. Des- 

barres, Eur. Polym. J. ,1979, 15, 993. 
[28] Gilpin, D. K.; and Geindoga, E., J. Chromatogr. 

Sci. ,1983, 21, 352. 
[29] Cosgrove, T.; Vincent, B.; Stuart, M. C.; Barnett, K.G.; 

and Sissons, D.S., Macromolecules, 1981, 14, 1018. 
[30] Daoud, M.; and Cotton, J. P., J. Physique, 1982, 43, 

531. 
[31] Witten, T. A.; and Pincus, P. A., Macromolecules, 1986, 

19, 2509. 
[32] Borukhov, I.; and Leibler, L., Phys. Rev. E, 2000, 62, 

R41. 
[33] Dan, N.; and Tirell, M., Macromolecules, 1992, 25,2890. 
[34] Toral, R.; and Chakrabarti, A., Phys. Rev. E, 1993, 47 



,4240. 
[35] Grest,G.S.; and Murat, M., Macromolecules 1993, 26, 

3108. 
[36] Wijmans, CM.; Leermakers, F.A.M.; and Fleer, G.J., 

Langmuir 1994, 10, 4514. 
[37] Gcrda, J. J.; Sintes, T.; and Toral, R., Macromolecules 

2003, 36, 1407. 
[38] Baungartner, A., in Applications of the Monte Carlo 

Method in Statistical Physics, 2"^* ed. Edited by K. 

Binder. Topics in Current Physics 36. Springer- Verlag, 

Berlin, 1987. 
[39] Allen, M.; and Tildesley, D., Computer Simulation of 

Liquids; Clarendon, Oxford, 1987. 
[40] de Gennes, P. G., Scaling concepts in polymer physics; 

Cornell University,Ithaca NY, 1979. 
[41] Edwards, S.F., Proc. Phys. Soc. ,1965,85, 613. 
[42] Garignano, M.A.; and Szleifer, I., J. Chem. Phys, 1995, 

102, 8662. 
[43] Flory, P. J. Principles of Polymer Chemistry; Cornell 

University Press, 1953. 
[44] Cosgrove, T.; Heath, T.; van Lent, B.; Leermak- 
ers, F.A.M.; and Scheutjens, J. M. H. M., Macro- 
molecules, 1987, 20, 1692. 
[45] Ruckenstein, E.; and Li, B., J.Chem.Phys 1997, 107, 3. 
[46] Lin, E.K.; and Cast, A. P., Macromolecules,1996,29, 390. 
[47] Evans, R.; and Napper, D. H., J. Colloid Interface Sci, 

1978, 63, 43. 




FIG. 1: Schematic representation of a spherical brush with 
an impenetrable core of radius Tc inside a spherical cavity of 
radius R. 
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FIG. 2: Comparison between the monomer density profile 
(f>{r) for free or uncompressed spherical polymer brushes ob- 
tained from our MC simulations (lines) and SCF calculations 
(symbols) for different values of the chain length A'^ and num- 
ber of grafted chains /. The results are compared to the 
predictions of Cariagno and Szleifer^J (CS, filled symbols) 
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FIG. 7: Log-log plot of the system pressure P vs. the 
monomer volume fraction v obtained from MC simulations for 
different values of the chain length A^ and number of grafted 
chains /. A solid line of slope 2.73 is included to guide the 
eye. 
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